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ABSTRACT 


Some  response  statistics  of  a  linear  system  under  impulsive 
forcing  are  here  examined  from  analytical  as  well  as  empirical 
viewpoints.  In  particular,  the  first-order  probability  density 
and  crest  and  rise  statistics  of  the  response  of  a  single-degree- 
of-freedom  system  under  Poisson  impulsive  noise  forcing  are 
considered.  The  failure  of  all  attempts  to  evaluate  the  Gram- 
Charlier  series  expansion  for  the  response  probability  density 
is  noted,  although  from  these  considerations  a  convenient  recur¬ 
sion  relation  between  the  cumulants  and  the  moments  of  a  random 
process  is  derived.  Simultaneous  assumptions  of  low  impulse 
frequency  and  high  system  Q  lead  to  analytical  determination  of 
the  response  crest  distributions,  which  are  also  verified  experi¬ 
mentally.  In  view  of  the  absence  of  applicable  theory,  the  rise 
statistics  are  examined  empirically  and  compared  with  the  crest 
statistics  obtained  under  like  conditions.  Finally,  damage  rate 
estimates  based  upon  the  aforementioned  statistics  are  obtained; 
typical  calculations  predict  the  damage  incurrence  for  a  linear 
svstem  under  Poisson  impulsive  noise  forcing  to  be  most  conserva¬ 
tive  when  crest  statistics  are  used  and  also  substantially  greater 
than  the  "worst  possible"  gaussian  case. 
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Section  I.  Introduction 


This  report  considers  some  statistical  aspects  of  the  response 
of  a  linear  single-degree-of-f reedom  system  to  non-gaussian  loading 
relevant  to  damage  and  reliability  considerations.  In  particular, 
the  response  first-order  probability  density,  the  crest  statistics, 
and  the  rise  (or  fall)  statistics  are  examined  empirically  and, 
whenever  possible,  analytically. 

The  non-gaussian  noise  sources  here  considered  in  detail  are 
those  which  may  be  closely  characterized  by  a  Poisson  impulsive 
noise  model.  Common  examples  of  impulsive  noise  are  gust  loading 
of  space  vehicles,  lightning-induced  static,  switching  noises  on 
transmission  lines,  and  ambient  noise  in  the  ocean.  These  can  be 
identified  either  by  a  highly-peaked  behavior  (compared  with  the 
normal  distribution)  at  the  mean  values  of  their  first-order  prob¬ 
ability  densities  or,  in  the  time  domain,  by  large  high-amplitude 
bursts  of  noise  against  low-level  backgrounds.  Unfortunately,  in 
an  analytical  sense,  system  response  to  impulsive  noise  forcing 
does  not  in  general  lend  itself  to  the  frequent  simplifications 
associated  with  gaussian  noise  forcing,  e.g.,  a  linear  operation  on 
a  gaussian  process  yields  another  gaussian  process.  However,  if  a 
reasonably  accurate  characterization  of  the  non-gaussian  process  is 
realized  by  using  a  Poisson  impulsive  model,  then  inherent  analytical 
difficulties  can  be  tolerated  and  hopefully,  with  better  understanding, 
overcome.  Consequently,  because  of  the  search  for  more  accurate 
models  of  certain  field  data,  an  increasing  amount  of  present-day 
interest  in  stochastic  analysis  of  engineering  systems  has  shifted 
toward  consideration  of  non-gaussian  processes,  such  as  impulsive 
noise . 

First,  we  examine  in  detail  theoretical  representations  of  the 
first-order  probability  density  of  the  response  of  the  single-degree- 
of-f  reedom  (1D0F)  system  to  impulsive  noise  forcing.  It  is  possible 
that  under  certain  conditions  the  system  excitation  becomes  a  train 
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of  non-overlapping  responses  (i.e.,  to  individual  impulses)  whereby 
the  response  problem  reduces  simply  to  that  of  determining  the 
probability  density  of  a  single  scaled  impulse  response  over  the 
inverse  of  the  average  frequency  of  occurrence  for  the  impulsive 
noise.  Rather  more  general  approaches  are  utilized  here;  in 
particular,  an  evaluation  of  the  Gram-Charlier  series  expansion  for 
the  response  probability  density  is  attempted.  A  general  recursion 
formula  relating  the  coefficients  of  two  series,  one  being  the 
logarithm  of  the  other,  is  found  and  then  is  used  to  evaluate  coeffi¬ 
cients  of  the  Gram-Charlier  series  in  terms  of  the  cumulants  of  the 
response  process.  This  relation  also  allows  one  to  express  the 
response  moments  in  terms  of  the  cumulants,  or  vice  versa,  in  a  com¬ 
pact  fashion.  Unfortunately,  it  was  found  that  usable  computations 
of  the  first-order  probability  density  for  comparison  with  experi¬ 
ment  cannot  be  performed  via  the  Gram-Charlier  series.  Computational 
difficulties  encountered  in  finding  the  probability  density  serve  not 
only  as  a  dire  warning  that  an  elegant  theoretical  result  may  not  be 
very  helpful  in  numerical  work,  but  also  as  a  possible  inducement  for 
other  investigators  to  develop  alternate  approaches  to  the  computa¬ 
tion  problem. 

Next,  the  crest  (level  crossing)  statistics  of  the  response  of 
the  1D0F  system  under  Poisson  impulsive  noise  loading  are  considered 
from  both  analytical  and  empirical  viewpoints.  Theoretical  computa¬ 
tions  are  here  obtained  for  high  Q  systems  driven  by  noise  with  low 
impulse  frequencies  (rates)  and  then  compared  with  laboratory  meas¬ 
urements.  In  addition,  some  empirical  response  probability  density 
curves  are  presented  which  establish  the  large  probabilities  expected 
for  the  occurrence  of  the  higher  stress  levels. 

In  the  third  section,  experimental  data  on  the  response  rise 
statistics  are  presented,  providing  a  measure  of  the  stress  differ¬ 
ences  between  response  minimum  and  the  next  maximum.  These  results 
are  then  compared  with  crest  statistics  for  which  approximations 
made  during  the  theoretical  development  are  valid  and  acceptable  for 
damage  predictions. 

Finally,  the  results  (analytical  and  empirical)  of  the  crest 
and  rise  studies  are  applied  to  sample  damage  rate  calculations 
using  the  Palmgren-Miner  cumulative  damage  law.  The  actual  numeri¬ 
cal  values  obtained  should  not  be  taken  too  literally  since  con¬ 
troversy  still  surrounds  fatigue  theory  for  random  loading.  Rather, 
it  is  the  trends  exhibited  in  these  calculations  by  varying  system 
Q  and  impulse  frequency  v  which  constitute  the  summum  bonum  of  this 
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section.  It  is  thought  that  the  tabulated  damage  rates  represent 
lower  bounds  on  the  more  precise  values;  thus,  one  infers  that  a 
gaussian  process  assumption  generally  represents  at  least  an  order 
of  magnitude  error  in  the  estimation  of  expected  system  lifetime 
under  impulsive  noise  forcing. 

Section  II .  Response  Probability  Density  for  Impulsive  Forcing 

Knowledge  of  the  first-order  response  probability  density  func¬ 
tion  of  a  1D0F  system  excited  by  non-gaussian  noise  is  of  interest 
in  nearly  all  predictions,  regardless  of  ultimate  applications.  It 
is  thus  logical  to  consider  this  response  function  first,  particularly 
since  portions  of  this  report  are  concerned  with  damage  predictions 
where  information  about  occurrence  of  large  stress  levels  is  most 
relevant.  After  establishing  certain  background  information,  the 
theoretical  investigation  examines  the  problem  of  obtaining  relation¬ 
ships  between  the  coefficients  of  two  power  series,  one  the  logarithm 
of  the  other,  resulting  in  a  simple  recursion  relation  between  the 
moments  and  the  cumulants  of  any  stochastic  process.  Next,  the 
Gram-Charlier  series  for  the  first-order  response  probability  density 
of  a  linear  system  is  examined  with  the  actual  numerical  calculation 
of  the  series  being  the  res  principalis.  Finally,  several  futile 
attempts  at  evaluating  the  response  probability  density  by  using 
numerical  integration  techniques  are  presented. 

A.  Background 

The  most  promising  theoretical  noise  model  for  engineering 
application  of  impulsive  noise  statistics  incorporates  a  Poisson- 
process  approximation  of  the  noise  source2-4  involving  sums  of 
impulse  (delta)  functions.  Although  from  a  strictly  mathematical 
viewpoint  the  Poisson  model  may  not  be  particularly  desirable,  to 
the  engineer  it  provides  a  reasonably  straightforward  and  intuitively 
appealing  attack  on  system  response.  Such  a  model  for  impulsive 
noise  enables  one  to  employ  the  convenient  properties  of  delta  func¬ 
tions  in  various  integral  transformations.  The  delta  function 
formally  may  be  defined  by  the  pair  of  relations5 


8(t)  =0,  t  ^  0 


(1) 


and 
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t  +T 

O 

g(t)8(t-t  )  dt  =  g  (t  ),  t  >  0  (2) 

o  o 

t  -T 

o 

for  any  function  g(t)  continuous  at  t  .  Of  course,  theoretical 
results  obtained  via  this  method  must  never  be  accepted  without 
question;  witness,  for  example,  some  of  the  paradoxes  encountered 
in  communications  theory  or  quantum  mechanics^  through  indiscreet 
use  of  the  delta  function.  However,  the  pathological  properties 
of  the  delta  function  can  be  kept  under  control,  if  one  keeps  in 
mind  the  meaning  of  the  limiting  processes  involved  in  its  defini¬ 
tion.^'®  If  proper  care  is  exercised,  one  obtains  theoretical 
results  which  are  physically  meaningful  and  which  approximate 
empirical  results  quite  well. 

The  Poisson  model  of  impulsive  noise  is  realized  as  a  train 
of  such  delta  functions,  namely, 


n(t)  =  £  a±  8  (t-t±)  .  (3) 

i 

Here,  t.  is  the  occurrence  time  of  the  ith  impulse  while  the  a. ’s 
are  the  impulse  strengths,  the  "areas"  under  each  delta  function. 
Now  in  order  to  obtain  tractable  results,  the  t. ’s  and  the  a.  *s 
must  conform  to  certain  statistical  requirements.  For  the  Poisson 
model  the  requirements  are  as  follows.  First,  the  occurrence 
times  are  assumed  to  be  independent  and  identically  distributed, 
having  a  Poisson  distribution. 


PT(n) 


(v.T)n 

n  J 


— VT 

e 


(4) 


for  the  probability  of  n  impulses  occurring  in  an  observation  time 
T  with  v  being  the  average  impulse  frequency.  Second,  the  model 
requires  the  a^ 1 s  to  be  independent  random  variables,  identically 
distributed  with  probability  density  p(a).  If  the  above  require¬ 
ments  for  t-  and  a-  are  fulfilled,  then  v  and  p(a)  characterize 
the  Poisson  noise  model.  ' 3 
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In  this  report  we  are  studying  the  response  of  a  single-degree- 
of-freedom  (1D0F) ,  strictly  a  "single  mode",  system.  If  the  forcing 
function  is  known  explicitly  as  a  function  of  time,  then  the  response 
of  the  1D0F  system  (and  hence  all  statistical  properties  of  it)  can 
be  obtained  from  solutions  to  the  (scaled)  standard  form  of  a  linear 
second-order  differential  equation, 

+  2 £o>  +  m  2  x(t>  =  f(t)  ,  (5) 

dt2  °  dx 

where 

f(t)  =  the  (scaled)  stochastic  forcing  function, 

x (t )  =  the  system  stochastic  response, 

and  £  and  |ge  the  damping  ratio  and  undamped  resonant  frequency, 
respectively.  Lacking  such  complete  description,  certain  sta¬ 
tistical  properties  x(t)  may  be  obtained  by  using  the  formal 
solution  to  Eq.  (5) ;  such  information  readily  yields  various 
response  probability  densities  when  f(t)  is  gaussian.  However,  for 
impulsive  noise  forcing  one  is  not  so  fortunate. 


Since  the  impulse  responsj2 function  h(t)  determines  the  system 
response  to  arbitrary  forcing,  one  ip^ght  envision  its  usefulness 
here.  Now,  for  an  underdamped  system 


h  (t) 


t  >  0 
t  <  0 


(6) 


where  cc  and  |3  are. 


and 


respectively, 

0) 

a  =  -2-=  (jo  r 

2Q  o 


p  =  00 


77T 


<-V 


(7) 

(8) 


The  parameter  Q  is  the  system  quality  factor  while  A  is  simply  a 
gain  factor,  implicitly  dependent  on  the  impulse  strength. 


Thus,  the 
system  response 
V,  p  (a)  ,  coq,  and 


performance  of  both  the  excitation  and  the  dynamic 
are  characterized  here  by  four  parameters,  namely 

Q. 
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B.  A  Recursion  Relation 

Let  us  now  digress  and  examine  the  problem  of  relating  the 
coefficients  of  two  pov/er  series,  one  the  logaritlmof  the  other. 
Consideration  of  this  seemingly  abstract  mathematical  problem  may 
at  first  seem  unwarranted,  but  its  solution,  when  applied  to  the 
problem  of  relating  the  nth  moment  to  the  nth  cumulant  of  a  stoch¬ 
astic  process  will  prove  most  useful.  Thus,  consider  the  equation 


g(x)  =  eJ 


(x) 


(9) 


where  g (x)  and  f (x)  have  the  power  series  representations 


oo 

(*>  -  Z9(n)(0)  £  , 


0 

oo 


f(x)  =  2  f(n)  (°)  77  . 

0 


n 

c 

n' 


(10) 


(ID 


The  problem  is  to  relate  the  coefficients  (o)  and  g^  (0)  . 

Equation  (9)  can  be  differentiated,  yielding 


g(1)  (x)  =  f  (1)  (x)ef  (x) 

-  f(1)  (x)g(o)  (x) 


(12) 


then,  the  successive  derivatives  are  given  by 


g(2)  (x) 
g(3)  (x) 
g(4)(x) 


f(2)  (x)  g  (0)  (x)+f(1)  (x)g(1)  (x) 

f  (3)  (x)g(0)  (x)+2f  (2)  (x)g(1)  (x)+f  (1)  (x)g(2)  (x) 


(13) 


f(4)  (x)g(0)  (x)+3f(3)  (x)g(1)  (x)+3f(2)  (x)g(2)  (x)+f(1)  (x)g(3)  (x) 


g(l)  (x)  =  ^  (j"3)  f(1' (x)g'— '  (x) 


.-k+1)  .  ,  (k-1) 


(14) 


k=l 
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(15) 


where 


is  the  usual  binomial  coefficient 


i-l\  =  (i-1)  : 

k-l)  (k-1) ! (i-k) ! 


If  Eq.  (14)  is  evaluated  for  zero  argument,  then 

g(l)(0)  =  £  (£li)  fU"k+1)  (0)g(k"1)  (0)  ,  (16) 

k=l 

which  is  exactly  the  desired  relationship  between  the  coefficients 
of  the  two  power  series  for  g (x)  and  f (x) . 


Now,  one  can  consider  the  problem  of  relating  the  nth  moment. 


m  ,  to  the  nth  cumulant,  X  .  Specifically,  one  has 
n  n  J 

Zo°  .  n 

(iu) 

X  3 — , 

=  e  1  n  ni 


V  m 

Li  n  n  I 


(17) 


0 


Equation  (16)  is  directly  applicable  so  that 


mn  ^  (k-l)  Xn-k+l  ""k-l  ' 

k=l 

and  noting  that  m  =  1 , 

^n  mn  X  (k-l)  ^n-k+1  ^V-l 
k=2 


(18) 


(19) 


which  completes  the  solution!  The  general  results  of  this  section 
will  be  applied  in  the  following  section  where  the  Gram-Charlier 
expansion  will  be  considered. 


C.  A  Series  Expression  for  the  Probability  Density 

In  this  section,  a  formal  theoretical  analysis  of  the  response 
probability  density  for  a  lDOF  system  under  Poisson  loading  is  pre¬ 
sented.  For  the  most  part,  this  work  is  a  recapitulation  of  currently 
available  theory  .2 However,  the  recursion  relation  developed  in 
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the  preceding  section  is  used  to  assist  in  the  computation  of  the 
coefficients  of  the  series  expression  obtained  for  the  output  prob¬ 
ability  density.^  In  addition,  evaluation  techniques  for  deter¬ 
mining  the  response  transform  of  the  probability  density  are 
considered. 

The  inverse  transformation  of  the  characteristic  function, 
expressed  in  exponential  form,  is  formally  given  by 


p(y) 


(20) 


where  (by  definition)  the  A.n's  are  the  cumulants  of  the  process  and 
y  is  the  response  variable  of  the  1D0F  system.  It  is  convenient  to 
convert  y  into  a  normalized  zero-mean  unity-variance  variable  T). 
Then,  using  the  identity 

/  •  \  n 


n  n .' 


cd  . .  .  n 
=  V  a 

L  n  n .' 


(21) 


0 

the  probability  density  of  T)  may  be  shown  to  be  given  by 

1 


P(T1)  = 


2 

-T]  /2 


/  2tt' 


oo  a  H  (r|) 

Y  —0.  n 

L  n  n ' 


(22) 


0 


where  is  the  nth  Hermite  polynomial, 


H  (x)  =  e 
n 


+x2/2 


r  d_i 

L- 


dnn  -x  /2 
e 


(23) 


Now  the  Xn's  in  Eq.  (21)  have  been  shown  to  be  given  for  Poisson 
loading  by  an  extension  of  Campbell's  Theorem,2  i.e.. 


X 

n 


(24) 
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For  a  linear  1D0F  system,  with  h(t)  given  by  Eq.  (6),  one  obtains. 


X 

n 


<an) 

<a2>n/2 


< 


n-2 


0 


n  odd  , 


where  A  has  been  chosen  (see  Eq.  (46))  so  that  the  variance,  X^# 
is  unity.  Thus,  with  a,  fJ ,  v,  and  p(a)  known,  the  response  prob¬ 
ability  density  of  the  1D0F  system  under  impulsive  loading  has  been 
calculated  in  a  formal  way.  (However,  insofar  as  actual  numerical 
computations  are  concerned,  the  utilization  Eq.  (22)  leads  to 
apparently  insurmountable  difficulties  as  will  be  explained.)  The 
cumulants  may  be  calculated  straightforwardly,  and  the  evaluation 
of  the  series  coefficients,  a  ,  can  be  performed  by  utilizing 
the  considerations  of  the  preceding  section. 

Thus,  applying  Ej.  (16)  to  Eq.  (21),  a  recursion  relation  may 
be  established  for  the  nth  coefficient  of  the  series  expansion 
Eq.  (22) .  In  particular,  one  obtains 

a  =1 
o 

al  =  a2  =  0  , 
n-2 

a  =  X  +  y 

n  n  L, 

k=4 

where  the  summation  is  collapsed  until  n  =  6.  Thus,  Eq.  (26)  along 
with  Eqs .  (22)  and  (25)  permit  analytical  evaluation  of  the  response 

probability  density  under  impulsive  loading  by  digital  computer 
methods.  Hence  it  would  seem  that  the  general  problem  is  for  all 
intents  and  purposes  solved.  Unfortunately,  such  is  not  the  case 
here,  as  will  be  discussed  shortly. 

One  desires  to  carefully  select  the  impulse  strength  distribu¬ 
tions  p (a)  so  as  to  closely  model  field  problems.  Thus,  if  the 


(k-l)  Xn-k+l  ^-1 


n  >  3 


(26) 
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noise  source  produces  constant  amplitude  bursts  of  symmetric  noise, 
a  logical  choice  for  the^mpulse  distribution  of  the  model  would 
be  (normalization  apart) 

p  (a)  =  j  S(a+1)  +  j  8  (a-1)  ,  (27) 


hereafter  referred  to  as  the  fixed  symmetrical  impulse  strength 
(or  more  compactly,  FSIS)  case.  On  the  other  hand,  if  the  noise 
source  produces  pulses  of  random  strengths,  then  a  continuous  dis¬ 
tribution  of  the  form 


P  (a) 


_ 1_ 

/2irl 


(28) 


might  be  employed,  abbreviated  as  the  GDIS  (gaussian  distributed 
impulse  strength)  case.  These  two  cases  will  be  considered  in 
detail.  anticipating  early  success,  computations  of  the  response 
probability  density  were  made  under  the  more  general  GDIS  assump¬ 
tion.  The  infinite  series  in  Eq.  (22)  was  summed  on  a  CDC  1604 
digital  computer  at  the  Numerical  Analysis  Center,  University  of 
Minnesota,  for  n  as  large  as  300  and  for  a  range  of  veilues  of  sj  . 

Careful  examination  of  the  partial  sums,  so  obtained,  revealed 
the  series  to  be  slowly  diverging,  presumably  implying  that  the 
partial  sums  were  not  extended  sufficiently  high  in  n.  However, 
since  the  computational  limits  of  the  1604  were  already  reached,  n 
could  not  be  increased.  Such  problems  eventually  led  to  considera¬ 
tions  of  the  FSIS  case  mentioned  above.  One  notes  that  the  moments 
of  the  gaussian  distribution.  Eg.  (28),  increase  with  n.-*-^  However, 
for  the  FSIS  distribution  of  Eq.  (27) ,  all  moments  are  unity,  im¬ 
plying  that  the  cumulants  for  the  FSIS  case  will  be  correspondingly 
smaller  than  those  produced  by  the  GDIS  condition.  Therefore, 
faster  convergence  of  the  probability  density  series  is  expected. 

Although  the  partial  sums  obtained  on  the  1604  under  the  FSIS 
assumption  were  comparatively  "well-behaved"',  it  was  only  for  large 
v  when  the  resultant  distribution  varied  only  a  few  percent  from 
the  gaussian  that  rapid  convergence  could  be  obtained.  Since  the 
gaussian  limit  for  high  v  is  a  well  known  result^ even  for  more 
general  Poisson  noise  models,  these  "almost-gaussian"  cases  pro¬ 
vided  no  new  insight  into  the  numerical  computation  problem.  It 
was  concluded  that  obtaining  probability  densities  of  the  more 
non-gaussian  cases,  which  includes  many  of  more  practical  interest, 
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were  inaccessible  by  the  Gram-Charlier  series  approach.  Thus, 
other  approaches  must  be  developed  for  computing  the  first-order 
response  probability  density  of  the  1DOF  system  under  low  density 
impulsive  noise  loading. 


Two  additional  methods  of  analysis  were  tried,  both  still 
under  the  FSIS  hypothesis.  In  an  alternate  approach,  one  uses  the 
result  that  any  symmetric  (i.e.,  even  function)  probability  density 
may  be  formally  written  as^ 


P(?7) 


/-j  2n  (2n)  J 
1 


-iuT7 
e  ' 


du 


(29) 


Noting  that  the  series  is  an  alternating  one  and  that  the  nth 
coefficient  of  u  will  eventually  approach  zero  in  the  FSIS  case, 
one  can  presumably  sum  this  series  to  any  given  accuracy. 17 
Unfortunately,  for  the  interesting  non-gaussian  cases,  the 
characteristics  function  for  p  (77)  ,  is  expected  to  remain  large 
(compared  to  the  gaussian  characteristic  function)  for  large 
values  of  the  argument  u.18  In  order  to  perform  the  inverse  trans¬ 
formation,  the  series  in  the  exponent  must  necessarily  be  summed 
for  values  of  u  »  5.  For  such  large  values  of  the  argument, 
however,  the  terms  of  the  series  initially  increase  with  n;  in 
fact,  sample  calculations  show^that  the  meignitude  of  the  largest 
term  may  be  of  the  order  of  10  o^more.  Since  the  final  sum  is 
expected  to  be  of  the  order  of  10  or  less,  it  is  rather  obvious 
that  many  significant  figures  must  be  carried  with  each  calcula¬ 
tion  in  order  to  insure  accuracy.  Inherent  computer  limitations 
provide  an  impasse  to  utilization  of  this  approach. 

Still  a  third  method  was  used  in  an  attempt  to  obtain  the 
first-order  response  probability  densities  analytically.  Noting 
that  the  general  representation  for  the  probability  density  of  a 
Poisson  pulse  process  is 
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Under  FSIS  conditions,  the  above  reduces  to 


oo 


1 


CO 


\ 


CD 


dt  [cos  [uh  (t )]  -  l] 


-iuij 


(31) 


du  e 


e 


-oo 


In  principle,  the  exponent  in  the  characteristic  function  may  be 
computed  by  direct  numerical  evaluation  of  the  t— integral •  Exam¬ 
ination  of  the  t-integrand,  however,  shows  a  rapidly  oscillating 
function  for  large  values  of  u.  Altnough  rather  sophisticated 
technicrues  were  employed  to  evaluate  this  t— integral ,  it  soon  be 
came  evident  that  an  accurate  value  of  the  characteristic  function 
could  not  be  obtained  without  consuming  considerable  computer 
time.  If  nothing  else,  one  is  led  to  an  appreciation  of  differ¬ 
ences  between  the  "theoretical"  and  the  "practical"  solution  of  a 
statistical  problem. 

With  these  three  different  numerical  methods  leading  to  abject 
failure  in  each  attempt,  it  was  concluded  that  alternate  representa¬ 
tions  for  the  probability  density  obtained  by  using  Campbell's 
Theorem  are  useful  only  for  the  closely  gaussian  cases;  for  the 
more  interesting  non-gaussian  cases,  their  mysteries  are  inaccessible. 
It  is  believed  that  the  three  aforementioned  procedures  were 
sufficiently  exhausted  to  warrant  the  above  generalization,  pending 
the  discovery  of  newer,  more  successful,  computational  techniques. 

Section  III, The  Crest  Statistics  for  Impulsive  Forcing 

In  initial  studies  of  random  processes,  the  level-crossing 
statistics  y8re  amon<?  the  first  to  be  examined  relative2to  ratigue 
prediction,  9  being  the  easiest  to  handle  theoretically  and 
experimentally.  Moreover,  in  certaig^cases, they  provide  an  upper 
bound  on  other  "peak"  distributions.  Thus,  in  work  with  impul¬ 
sive  loading  an  investigation  of  the  crest  statistics  of  linear 
system  response  is  also  of  interest  and  will  be  presented. 

The  particular  analytical  analysis  of  the  crest  statistics 
considered  here  may  be  classified  as  an  "integral  approach  to  the 
problem.  This  does  not  preclude  other  analyses;  for  example,  a 
series  approach  similar  to  that  used  for  the  first-order  probability 
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density  might  be  considered.  Suffice  it  to  say  that  the  inability 
to  sum  earlier  series  encourages  one  to  here  employ  a  "non-series" 
approach.  Thus,  to  assist  in  the  theoretical  development,  it  is 
expedient  to  restrict  the  discussion  to  1D0F  systems  with  high  Q 
(say  Q  >  5)  and  to  noise  sources  with  low  values  of  impulse  fre¬ 
quency  v  (say  v  <  100) .  These  two  assumptions  are  not  unrealistic 
or  overly  restrictive  and  certainly  have  considerable  importance 
in  response  studies  involving  electro-mechanical  systems. 

A.  Theory 

Let  us  consider  the  physical  significance  of  the  above  restric¬ 
tions.  The  high  Q  approximation  assures  many  cycles  of  the  system 
response  during  one  time-constant  of  the  system's  envelope  decay 
transient.  The  second,  that  of  a  low  v,  implies  that  on  an  average 
the  impulses  occur  infrequently  enough  so  as  to  allow  a  nearly  com¬ 
plete  decay  of  system  response  before  the  occurrence  of  the  next 
pulse.  Now  one  must  keep  in  mind  the  interdependence  of  these 
assumptions;  the  second  can  only  be  wisely  used  after  the  first  has 
been  selected.  For  if  the  conditions  are  applied  in  the  reverse 
order,  then  given  a  certain  impulse  density,  one  might  choose  Q  too 
large  causing  the  system  response  to  not  decay  sufficiently  before 
the  occurrence  of  the  next  pulse.  Under  the  above  hypothesis,  a 
typical  time-domain  response  of  the  1D0F  under  impulsive  loading 
would  be  that  shown  below. 
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Analytical  determination  of  the  crest  statistics  of  a  stoch¬ 
astic  process  yields  knowledge  of  the  number  of  crossing  of  a  level 
77  as  a  function  of  that  level.  For  the  problem  under  considera¬ 
tion  this  can  be  accomplished  by  first  noting  that  the  high  Q 
assumption  assures  that  whenever  the  envelope  of  the  response  of 
the  system  under  study  is  above  the  level  rj  ,  each  level  crossing 
will  occur  at  a  rate  essentially  determined  by  the  resonant  fre¬ 
quency  f  of  the  system.  Thus,  the  number  of  crossings  at  17  is 
the  product  of  f  and  the  percentage  of  time  the  envelope  spends 
above  7 7  .  But  the  latter  is  simply  the  probability  of  exceedance 
of  the  envelope.  Hence,  we  write 

00 

N+  (77)  =  fQ  x  j  P£(§)  ,  (32) 

where  N+ (77)  is  the  number  of  positive  slope  crossings  of  the  level 
77  ,  and  the  integral  is  the  probability  of  exceedance  of  the  same 

level  77  by  the  envelope.  For  convenience,  Eq.  (32)  ijiay  be  normal¬ 
ized  to  the  number  of  zero  crossings  N  (0)  .  Using  c  (77)  to  denote 
the  normalized  crest  distribution,  one  has 


c+  (77) 


00 


pe(S)  as 


(33) 


Thus,  the  problem  of  crest  statistics  now  rests  with  finding  the 
probability  density  of  the  envelope,  Pv  (77)  . 


This  may  be  accomplished  by  using  a  reasonably  straightforward 
reinterpretation  of  S.  0.  Rice’s  integral  representation  for  prob¬ 
ability  density,  Eq.  (30).  Let 


c+  (77) 


(iu)  e 


-1UX 


du 


-00 


(34) 


here,  M  (iu)  is  the  transform  (i.e.,  characteristic  function)  of 
p  (77).  ENext,  one  argues  that  if  h  (t)  denotes  the  "impulse  response 
for  the  envelope  of  the  response  variable,  the  characteristic  func¬ 
tion  M  (iu)  can  be  written  formally  as 
E 
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under  the  stated  high  Q,  low  v  assumptions.  (An  explanation  for 
the  absolute  value,  |a| ,  follows  later.)  Thus  the  integral  repre¬ 
sentation  of  the  crest  distribution  becomes 


Although  this  equation  is  not  readily  integrable,  it  is  of  a  form 
to  which  the  third  approach  discussed  earlier  may  be  applied. 

For  sufficiently  small  impulse  densities,  v,  one  may  use  the 
relation 


e  '  =  1  + 


Ipl  <<  1. 


(37) 


Thus,  Eq.  (36)  becomes 

oo  oo 

c%)=^/  dx/  du  e— 


-oo 


+  -  /  dx 


-  EF  '  dx 


oo 


CD 


-CD 


CD  OO 


du  e 


-iux 


dt 


iu  |a|  h  (t) 
da  e  p(a) 


'O 


-oo 

oo 


du  e  1UX  I  dt  /  da  p(a)  . 


(38) 


-oo 


The  lower  limit  of  the  t-integration  becomes  zero  because  of  the 
physical-realizability  condition. 


h^  (t)  =  0  for  t  <  0 


(39) 
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By  noting  that  the  delta  function  may  be  expressed  as 


8(x) 


,  -iux 
du  e 


-oo 
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(40) 


the  first  and  the  last  terms  of  Eq.  (38)  can  be  formally  arranged 
to  give 


Let  us  examine  this  formulation  closely  for  a  proper  interpretation. 


The  t-integration  of  the  first  term  of  Eq.  (41)  appears  to 
render  the  method  useless.  However,  if  one  restricts  rj  to  be 
strictly  greater  than  zero, 


Tj  >  0 


(42) 


and  then  performs  the  x-integration  first,  the  first  term  of 
E1#  (41)  will  vanish  since  the  x-integration  involves  only  the 

strictly  zero-portion  of  the  delta  function.  Although  this  procedure 
of  eliminating  the  indeterminacy  of  the  t-integral  may  be  unorthodox, 
let  us  proceed. 


Equation  (41)  now  gives  the  crest  distribution  as 


c+  {r)) 


.tr«  J4 

d  o'  0 


OO 


27 r 


dt  /  da  e 


-iu  [x  -  |a|  h  (t)j 


p(0.  (43) 


oo  **  0  ^oo 


If  one  performs  the  u-integration  and  next  conducts  the  x-integra¬ 
tion,  the  crest  distribution  becomes 


r  t 

'  (rj)  =  V  I  dt  < 

J 0  ^oo 


da  p  (a)  x 


1  r)  >  |a|  hE  (t) 
0  otherwise 


(44) 


which  suggests  closer  examination  of  h  (t) . 

£j 
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Recall  that  one  is  here  working  with  envelope  statistics. 
Now  for  a  1D0F  system  under  the  high  Q,  low  v  hypothesis,  the 
envelope  "impulse  response"  is 


h£  (t)  =  A  e  at  t  >  0 


(45) 


where  the  gain  factor  A  is  the  same  as  that  of  Eg.  (6)  (for  normal¬ 
ization  purposes) .  By  setting  the  variance  to  unity,  the  gain 
factor  becomes 


A 


(— ) 
a  VB  ' 


1/2 


+  2‘ 


,  2\ 

<a  > 


(46) 


This  implies  an  equivalence  of  the  upper  inequality  of  Eq.  (44)  to 
the  inequality 


t  < 


la  I  A 


(47) 


But  since  the  time  range  of  interest  is 


t  >  0 


(48) 


Eq.  (47)  also  implies  that 


a 


(49) 


Therefore,  using  these  limits  and  performing  the  t-integration ,  the 
crest  distribution  becomes 


oo 


da  p (a)  log  — 

da  p  (a)  log  ~~ 


(50) 


Before  attempting  evaluation,  let  us  first  clarify  the  signifi¬ 
cance  of  the  second  integral  on  the  right-hand  side  of  Eq.  (50) ,  a 
somewhat  subtle  point.  The  envelope  response  of  a  1DOF  system  to 
an  isolated  positive  unit— amplitude  impulse  is  a  decaying  exponential 
always  greater  than  zero.  In  a  like  manner,  the  response  of  the 
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envelope  to  a  negative  unit  impulse  would  be  a  decaying  exponential, 
but  also  positive !  (This  explains  the  absolute  value  signs  placed 

around  the  impulse  strength  preceding  the  "envelope  impulse  response" 
h  (t) .)  One  must  keep  in  mind  that  he  is  dealing  with  envelope  sta¬ 
tistics  and  thus  their  attendant  subtleties.  Hence,  every  pulse, 
positive  or  negative,  will  produce  a  positive  envelope.  It  must 
also  be  kept  in  mind  that  overlapping  "envelope  responses"  due  to 
individual  pulses  do  not  add  to  give  total  envelope;  hence,  the  low 
V  and  high  Q  assumptions  are  important.  Restricting  the  analysis 
to  the  case  where  the  driving  impulses  have  a  symmetric  distribu¬ 
tion,  one  writes 

P  (a)  =  p(-a)  .  (51) 

Here,  the  average  occurrence  of  rates  of  the  positive  and  negative 
pulses  with  the  same  strengths  will  be  the  same.  Hence,  Eq.  (50) 
can  be  rewritten  as 

oo 

c+  (77)  =  ~  j  da  p  (a)  log  —■  .  (52) 

°Va 

The  above  is  a  form  useful  for  studying  the  crest  statistics  of  a 
1D0F  system  under  both  FSIS  and  GDIS  impulsive  loading. 


As  an  extra  bonus,  one  can  obtain  the  envelope  probability 
density  quite  easily  with  the  aid  of  Eq.  (50) .  From  Eq.  (33) , 

V7'1  =  -  %  c+(i»)  •  l53) 

By  applying  Eq.  (53)  to  Eq.  (52),  one  obtains 


pE(’)) 


00 


da  p (a) 


(54) 


as  an  expression  for  the  probability  density  of  the  envelope  of 
the  response. 


In  concluding  this  theoretical  development  of  the  response 
crest  statistics  of  a  1D0F  system  to  impulsive  loading,  it  is 
instructive  to  review  in  more  detail  the  restrictions  imposed  dur¬ 
ing  the  derivation  and  their  empirical  significance.  The  approxi¬ 
mation  of  Eq.  (37),  i.e.,  that  of  a  low  v,  certainly  presents  no 
experimental  problem;  thus,  practically  any  desired  accuracy  of  the 
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estimate  may  be  achieved.  Secondly,  one  might  ask  to  what  degree 
must  T|  be  greater  than  zero.  As  will  soon  become  evident,  (c.f. 

Eq.  (55))  there  are  some  cases  for  which  the  crest  distribution 
given  by  Er.  (52)  becomes  infinite  (not  unity)  as  Tf-*-  0.  However, 
if  acceptable  results  are  obtained  only  for  the  larger  stress  levels, 
this  is  sufficient  for  some  applications,  e.g.,  fatigue  studies  in¬ 
volving  an  endurance  limit.  The  pathological  behavior  at  the  smaller 
response  levels  can  thus  be  tolerated,  if  not  entirely  ignored. 

B.  Empirical  Results 

For  experimental  verification,  the  impulse  strength  distribu¬ 
tion,  namely  FSIS  and  GDIS  cases  considered  earlier,  were  investigated. 
In  order  to  carry  out  these  studies,  it  was  first  necessary  to  design 
an  impulse  noise  generator  whose  essential  characteristics  parallel 
the  Poisson  model,  and  then  to  select  a  linear  1DOF  system  whose 
response  can  be  studied  accurately.  The  generator  is  considered  in 
detail  in  Appendix  I.  Its  impulse  frequency,  v,  can  be  continuously 
varied  from  values  as  low  as  one  per  second  up  to  one  thousand  per 
second,  while  the  impulse  strength  distribution,  p(a),  can  be  set 
equal  to  the  first-order  density  of  any  process  which  may  be  conven¬ 
iently  generated  in  the  laboratory.  The  physical  system  chosen  for 
study  was  an  existing  electromechanical  device,  specifically  a  mass- 
loaded  cantilever  beam  driven  by  a  magnetic  field  transducer  and 
designed  to  simulate  a  1DOF  system.  This  device  has  been  described 
in  detail  in  other  published  works;  ^lt22  here,  suffice  it  to  say 
that  its  resonant  frequency  is  124.7  cps  and  the  system  quality 
factor  Q  can  be  varied  continuously  between  values  of  5  and  62.  This 
impulsive  noise  source  and  linear  system  were  used  for  all  experi¬ 
mental  work  reported  herein. 

Specializing  the  theory  worked  out  above  to  the  impulse  strength 
distribution  satisfying  the  FSIS  condition,  the  analytical  crest 
distribution  becomes 
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The  uniformity  in  the  appearance  of  the  product  cd  /qv  suggests  that 
this  ratio  best  characterizes  system  performance.0  However,  for 
these  studies  is  fixed  (784  to  be  exact)  ,  which  leaves  the  product 
Qv  to  be  examined.  Thus,  with  respect  to  this  particular  (JO  , 

Eq.  (55)  predicts  that  different  cases  of  Q  and  v,  but  with°constant 
Qv  products,  lead  to  the  same  theoretical  value.  This  result  was 
observed  (within  experimental  scatter)  in  all  experimental  work 
(probability  density,  crests,  and  rises);  sample  crest  data  appear 
in  Fig.  1.  Thus,  in  data  plots  only  a  single  experimental  set  of 
points  will  be  presented  for  each  Qv  value  and  will  represent  the 
average  of  several  runs. 

On  the  basis  of  experimental  work,  one  observes  that  rather 
drastic  departures  from  gaussian  process  statistics  are  readily 
obtained.  This  is  evident  by  examination  of  the  response  probability 
densities  (measured  simultaneously  with  the  crest  distribution) 
shown  in  Fig.  2  for  three  Qv  product  values.  They  indicate  prob¬ 
ability  densities  of  the  order  of  a  thousand  times  that  expected 
for  gaussian  densities  at  the  smaller  Qv ' s  at  response  levels  5  times 
the  rms  level.  The  Qv  product  provides  a  convenient  index  of  "non- 
gaussian"  performance. 

Comparisons  of  theoretical  response  crest  statistics,  Eq.  (55), 
with  experimental  crest  distributions  appear  in  Figs.  3,  4,  and  5. 
Assumptions  made  in  the  theoretical  derivation  suggest  that  better 
agreement  between  theory  and  experiment  will  be  achieved  for  lower 
values  of  v#  or  in  effect,  smaller  Qv  products.  Indeed,  this  was 
observed.  The  smallest  product  (Qv  =  38)  produced  results  in  best 
agreement  with  theory,  and  even  for  Qv  =  76,  only  slight  differences 
are  noted.  For  larger  Qv  products,  say  several  hundred  or  so,  the 
theory  is  not  valid  as  evidenced  in  Fig.  5.  However,  in  such  cases, 
e.g..  Fig.  5,  the  impulsive  noise  response  is  quite  gaussian  and  utili 
zation  of  a  gaussian  distribution  is  very  realistic.  Thus,  it  appears 
for  crest  statistics  under  the  FSIS  condition,  that  one  is  justified 
on  empirical  grounds  in  using  Eq.  (55)  for  Qv  <  100  (Qv/00  <  .128) 

and  the  gaussian  crest  distribution  for  larger  Qv  products. 

During  the  course  of  the  above  investigation,  certain  trends  in 
the  system  response  were  observed.  First,  as  the  impulse  rate  was 
increased,  the  response  statistics  approached  those  of  a  gaussian 
random  process;  it  has  been  shown  analytically-1-4  that  as  V  — ►  oo  ,  the 
first-order  probability  density  approaches  the  gaussian  distribution. 
Intuitively,  one  would  also  expect  the  crest  statistics  to  do  likewise 

Second,  intuitive  argument  (which  may  be  described  as  a  "folk 
theorem")  suggests  that  the  statistics  of  any  " reasonably-random" 
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Fig.  1  Comparison  of  FSIS 

Crest  Data  for  Constant 
Qv  Product  Employing  Q 


00 

lO 

o 

u 

II 

« 

II 

ii 

40 

« 

a 

B 

IT J 

H 

00 

00 

CD 

« 

CD 

ro 

in 

N 

Oi 

N 

— 

ITJ 

n 

n 

1 1 

n 

» 

o 

o 

O 

(0 

> 

O 

»• 

r* 

r* 

T3 

j- 

j- 

j— ' 

C 

ITJ 

O. 

O. 

CL 

X 

X 

X 

LU 

UJ 

llJ 

□ 

<1 

o 

o 

iri 


O 


o 

rd 


o 

oo 


{L)+o  S0I1SI1V1S  !S3dO 


21 


NORMALIZED  RESPONSE  VARIABLE 
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Fig.  2  The  Response  Probability  Density 
under  the  FSIS  Condition 
Employing  Qv  as  a  Parameter 
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process  when  passed  through  a  sufficient  narrow  band  filter  (i.e., 
a  high  Q  -  1DOF  system)  approach  those  of  a  gaussian  process.  When 
the  data  was  examined  as  a  function  of  Q  (i.e.,  v  held  constant) 
such  effects  were  clearly  evident.  Thus,  for  either  v  or  Q  suffi¬ 
ciently  large  (say,  according  to  Qv  >  500  or  Qv/o>  >  .64)  the  normal 
distributions  may  justifiably  be  used. 

Proceeding  to  the  more  general,  GDIS  condition  in  an  attempt 
to  further  verify  the  credibility  of  Eq.  (52)  ,  it  was  found  that 
solutions  of  Eq.  (52)  with  p(a)  gaussian  can  be  obtained  only  by 
numerical  integration.  Examination  of  the  problem  again  shows  that 
it  is  the  product  Qv  (with  Q  »  1)  which  controls  the  response  dis¬ 
tributions.  Paralleling  the  preceding  study,  crest  statistics  and 
first-order  probability  densities  obtained  experimentally  were 
nearly  identical  for  identical  Qv  product  values.  Thus,  "average" 
data  points  are  again  presented. 

Since  the  GDIS  case  is  of  more  general  interest  than  the  FSIS 
condition,  extended  investigations  were  carried  out.  Values  of 
the  Qv  products  as  low  as  Qv  =  7.6  were  obtained.  For  this 
"extremely  non-gaussian  case",  normalized  stress  levels  as  high  as 
10  times  rms  were  observed,  far  exceeding  those  obtained  in  earlier 
work;  detail  is  provided  in  Fig.  6  where  the  measured  probability 
density  for  Qv  =  7.6  is  tabulated.  A  discussion  of  these  prob¬ 
ability  density  measurements  could  parallel  that  of  the  FSIS  case 
presented  earlier,  but  with  one  notable  exception.  Here,  the 
probability  of  occurrences  at  the  higher  response  levels  is  con¬ 
siderably  greater  by  virtue  of  the  fact  that  the  impulse  strengths 
can  be  arbitrarily  large.  Results  of  both  numerical  integration 
of  Eq.  (52)  and  also  from  experimental  work  yield  crest  statistics 
lying  far  above  the  corresponding  gaussian  distribution.  Figures 
7,  8,  and  9  present  typical  examples. 

It  may  be  noted  that,  for  Qv  =  380,  the  experimental  and  theor¬ 
etical  crest  statistics  are  still  essentially  compatible.  Compare 
this  with  the  performance  under  the  FSIS  condition  for  the  same  Qv 
value.  One  therefore  concludes  that  the  theoretical  predictions  of 
Eq.  (52)  are  acceptable  under  the  GDIS  condition  even  for  high  Qv 
values,  and  this  in  spite  of  the  fact  that  the  GDIS  condition  is  a 
more  complex  case.  Hence,  for  values  of  Qv  up  to  say  400  (again, 

Qv/a30  <  *5)'  Ec2*  (52)  appears  valid;  for  larger  Qv’s,  one  may  use 
the  gaussian  distribution  with  but  little  observable  loss  in  accuracy. 
Similar  trends  were  observed  in  earlier  impulsive  noise  response 
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.  .  22 
studies,  utilizing  a  dead-zone-device  impulsive  noise  generator, 

though  such  data  do  not  seem  to  be  readily  interpretable  quantita¬ 
tively  in  terms  of  the  Poisson  noise  model. 

Finally,  a  comparison  of  the  crest  statistics  for  the  FSIS 
and  GDIS  conditions  for  identical  Qv  (=  38)  is  made  in  Fig.  10. 

Note  that  the  FSIS  prediction  is  similar  to  the  GDIS  prediction  at 
the  smaller  levels,  but  departs  rapidly  (actually  going  to  zero  at 
Tj  =  6.42)  from  the  GDIS  condition  which  continues  on  over  the  higher 
stress  levels.  Consequently,  any  damage  calculations  based  upon 
the  GDIS  model  will  certainly  predict  greater  damage  rates  than 
the  FSIS  model  at  the  same  rms  stress  level. 

Section  iv.  The  Rise  Statistics  for  Impulsive  Forcing 

This  section  will  consider  response  rise  statistics  of  a  1DOF 
system  subject  to  Poisson  impulsive  noise  loading.  Theoretical 
predictions  of  the  rise  statistics  (i.e.,  the  statistics  of  the 
diffg^ence  between  a  maximum  value  and  its  preceding  minimum) in 
work  to  date  have  been  limited  to  gaussian  processes.  Examina¬ 
tion  of  available  results  reveals  that  at  best  rather  restricted, 
approximate  solutions  to  the  rise  problem  have  been  obtained;  in 
particular,  only  gaussian  processes  with  rectangular  power  spectra 
have  been  considered.  Thus,  for  impulsive  noise  excitation,  present 
state-of-the-art  confines  one  t^4  experimental investigations , 
paralleling  the  work  of  Leybold  and  Fuller  for  gaussian  processes. 
Hopefully,  through  these  efforts,  a  better  understanding  of  the 
problem  may  be  achieved  v/hich  may,  in  turn,  provide  guidance  for 
further  analytical  developments. 


Since  the  £;Ls^g  are  related  to  the  "stress  reversals"  of 
fatigue  theory,  '  it  is  of  considerable  interest  to  use  both 
the  impulsive  noise  generator  and  the  1DOF  system  whose  crest  sta¬ 
tistics  were  discussed  in  the  preceding  section  so  that  crest  and 
rise  statistics  may  be  compared.  Of  considerable  importance  is  the 
development  of  a  rise-and-fall  distribution  analyzer  which  can 
operate  over  the  amplitude  and  frequency  rartges  of  interest.  Such 
an  electronic  installment  would  supersede  optical  methods  employ^ 
in  earlier  Judies  or  the  computer  processing  used  by  Leybold 
and  Fuller.  Appendix  II  describes  a  counting  device  recently 
developed  in  this  laboratory.  Operating  on  a  di f ferentiator-gated- 
integrator  principle,  it  processes  band-limited  signals  whose 
ratios  at  upper- to- lower  cutoff  frequencies  are  about  20:1  and 
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measures  rises  (their  amplitudes  being  normalized  with  respect  to 
twice  rms)  in  the  range  from  0.2  to  6.0.  In  fact,  if  small  rises 
situated  on  top  of  large  response  peaks  are  considered  unimportant, 
then  even  larger  rises  can  be  measured  without  additional  modifica¬ 
tion. 


Measurements  were  conducted  under  the  GDIS  condition,  this  case 
being  more  applicable  to  actual  field  data.  Many  values  of  Q  and  v 
were  studied  with  particular  attention  focusing  upon  comparison  of 
cases  with  like  Qv  products.  As  was  the  case  for  both  the  crest 
statistics  and  the  first-order  probability  density,  situations  with 
like  Qv  products  produced  nearly  identical  rise  statistics.  Again, 
without  loss  of  information,  data  shown  in  Fig.  11  represent  averages 
of  data  actually  recorded. 

In  earlier  studies  of  rise  statistics  of  several  filtered 
gaussian  processes,  it  was  found  that  rise  exceedance  probability 
distributions  were  bounded  above  by  the  integrated  Rayleigh  (the 
crest  distribution)  and  were  actually  bounded  by  positive  maxima 
statistics  (exceedance  probability  of  posi^ve  maxima)  for  processes 
with  identical  rms  bandwidth  parameter  e  .  Rise  distributions 
here  reported  for  Poisson  impulsive  noise  loading  of  a  1DOF  system 
behave  as  did  the  response  first-order  probability  and  crest  sta¬ 
tistics,  i.e.,  the  rise  statistics  fell  below  the  integrated  Rayleigh 
over  the  lower  stresses  and  lie  far  above  for  larger  values. 
Intuitively,  then,  one  expects  greater  damage  to  be  predicted  by 
using  rise  statistics  with  a  system  driven  by  impulsive  noise  than 
for  the  worst  possible  gaussian  case  at  the  same  rms  response  level. 
This  agrees  qualitatively  with  earlier  estimates  of  relative  damage 
based  on  crest  statistics;  both  of  these  estimates  will  be  more 
thoroughly  investigated  in  the  following  section.  More  important, 
we  note  that  for  Poisson  impulsive  noise,  the  crest  distribution 
provides  an  upper  bound  on  the  rise  exceedance  probabilities  for 
the  cases  studied;  see,  for  example.  Fig.  12. 

Figure  11  also  seems  to  show  that  the  rise  statistics  tend  to 
approach  an  integrated  Rayleigh  distribution  as  Qv- oo  .  However, 
this  limiting  result  is  not  true  in  general.  For  example,  if  one 
fixes  Q  very  small  and  then  lets  v~ oo  ,  the  system  response  will 
be  a  wideband  gaussian  process  whose  rise  statistics  will  fall  below 
the  integrated  Rayleigh.  Compare  this  with  v  being  held  small  and 
letting  oo  .  In  this  latter  case,  the  gaussian  process  is 

again  approached,  but  is  now  narrow-band  with  Rayleigh-distributed 
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Fig.  11  The  Response  Rise  Statistics 
under  the  GDIS  Condition 
Employing  Qv  as  a  Parameter. 
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23,27 


rises.  Thus,  the  respective  rise  statistics  are  different,  and 
the  two  limit  operations  would  approach  different  functional  forms 
appropriate  to  gaussian  processes  with  the  given  power  spectra, 

Such  differences  were  not  particularly  noticeable  empirically  since 
the  lowest  bound  on  Q  for  the  studies  was  7.6,  implying  that  the 
1D0F  system  still  behaved  essentially  as  a  narrow-band  filter. 

Hence,  for  our  limited  range  of  system  parameters,  the  rise  sta¬ 
tistics  would  indeed  limit  to  the  integrated  Rayleigh  as  Qv  ►  oo . 
For  low  values  of  Qv  very  radical  departures  from  the  Rayleigh  pre¬ 
dictions  are  empirically  obtained  at  the  "tail"  of  the  distribution, 
indicating  how  inaccurate  the  assumption  of  a  gaussian  response 
would  be  in  such  situations. 


Section  V.  Damage  Estimates  for  a  Linear  System 

The  motivation  for  studying  "peak"  statistics  in  connection 
with  aerospace  system  vibration  is  a  need  to  predict  system  damage 
due  to  wearout  by  acoustical  fatigue.  It  has  been  found  that  pre¬ 
diction  of  system  lifetime  is  a  rather  complicated  matter;  various 
theories  have  been  evolved^® but  the  only  tractable  and  widely- 
used  method  for  obtaining  damage  estimates  is  the  so-called 
Palmgren-Miner  law  of  linear  damage  accumulation.  This  law  will 
be  used  here  for  comparisons  between  damage  rates  for  the  linear 
harmonic  oscillator  under  impulsive  noise  forcing  and  for  gaussian 
forcing  with  the  same  rms  response. 

Thus,  for  comparison  purposes,  damage  estimates  will  be  based 
upon  a  linear  accumulation  law  with  damage  increments  independent 
of  the  order  in  which  the  "peak"  values  occurred.  Moreover,  it 
will  be  assumed  that  some  empirical  S-N  curve  is  available  and  as 
such  may  be  empirically  fitted  to  a  power  law  of  the  form 

S  ao 

N  ^  =  (— )  (56) 

N  v  S 

o 

where  aQ  usually  lies  between  5  and  15.  In  particular,  inter¬ 
dependence  of  stress  cycles  can  be  accounted  for  by  "rotating  the 
S-N  curves",  i.e.,  by  choosing  a  different  aQ  from  that  measured 
for  sinusoidal  stress.  Finally,  damage  calculations  will  be  .based 
upon  "peak  statistics"  as  given  by  the  crest  and  rise  statistics 
of  the  response  of  the  system  to  impulsive  loading.  The  predicted 
damage  per  second  using  each  of  these  definitions  of  "peak"  will 
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be  compared  with  the  damage  predicted  by  the  R  yleigh  distribution 
for  the  same  rms  response. 


First,  the  failure  rates  predicted  by  using  the  crest  statistics 
of  the  response  of  a  linear  single-degree-of-f reedom  system  leave 
out  double  maxima  and  other  minor  fluctuations  in  the  response. 

As  mentioned,  the  damage  rates  can  be  compared  to  those  given  by  a 
Rayleigh  distribution,  i.e.,  the  worst  possible  gaussian  cfise.  By 
denoting  these  failure  rates  as  <1/N)  and  <1/N>  ,  one  has 
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which  after  application  of  Eq.  (54)  becomes 
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Setting  =  5 ,  a  lower  bound,  damage  r.ites  for  several  Qv  products 
are  given  in  the  following  table  for  both  the  FSIS  and  GDIS  distribu¬ 
tions.  In  addition.  Table  I  has  the  damage  rates  of  E  t .  (57)  for 
the  same  obtained  by  numerical  integration  of  experimental  data. 


TABLE  I 


Predicted  Relative  Damage  Rates  for  a  1D0F  System  under  Impulsive 

Noise  Loading. 


Crest  Stat. 

Rise  Stat. 

Qv 

FSIS 

Cond . 
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Exp 
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43.2 

41.9 

76. 
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Of  course,  one  can  also  use  the  more  appropriate  rise  statistics, 
which  represent  the  "stress  reversal"  Statistics ,26  to  calculate 
damage  rates.  Since  no  useable  analytical  approach  to  the  rise 
statistics  has  been  developed  to  date, 23, 26  only  numerical  integra¬ 
tion  of  the  measured  rise  distributions  is  reported  here.  Statisti¬ 
cally  speaking,  data  points  at  the  higher  stress  levels  are  based 
on  small  samples  and  hence  are  not  too  accurate.  Since  the  power 
law,  Eq.  (56) ,  will  tend  to  emphasize  these  inaccuracies  at  the 
higher  response  levels,  failure  rates  obtained  via  this  method  will 
indubitably  be  inaccurate  in  actual  value,  but  will  presumably 
indicate  lower  bounds.  Thus,  for  several  Qv  products,  damage  rates 
based  on  the  fifth-power  law  are  shown  in  the  right-hand  column  of 
Table  I. 

It  is  noted  that  all  calculations  in  Table  I  show  decreasing 
damage  rates  for  a  1DOF  system  under  impulsive  noise  loading  as  Qv 
increases.  The  importance  of  knowledge  of  the  Qv  product  in  de¬ 
scribing  response  statistics  has  been  stressed  in  earlier  sections  and 
is  reaffirmed  by  these  damage  rate  calculations.  Moreover,  the  rates 
based  upon  the  theoretical  crest  statistics  (GDIS)  are  more  conserva¬ 
tive  (except  for  Qv  =  380  where  they  are  nearly  identical) .  (Compari¬ 
son  between  the  GDIS  experimentally-derived  damage  rates  and  those 
obtained  from  the  rise  statistics,  however,  show  one  smaller  estimate 
of  damage  for  the  former  probably  due  to  the  truncation  error  in  damage 
computation.)  Thus,  it  would  appear  that  fatigue  calculations  could 
be  based  on  theoretical  crest  statistics  in  order  to  obtain  a  con¬ 
servative  estimate  of  system  damage,  in  the  absence  of  rise  distribu¬ 
tions  (which  are  more  difficult  to  measure).  Nevertheless,  even  the 
damage  rate  predictions  based  empirically  on  the  rise  distributions 
generally  predict  an  order  of  magnitude  increase  over  the  appropriate 
Rayleigh  case,  excepting  say  when  Qv  >  380  (Qv/^q  >  0.5) . 

Comparison  between  the  FSIS  and  GDIS  columns  in  Table  I  also 
shows  the  importance  of  higher  stress  levels.  As  discussed  pre¬ 
viously,  crest  distributions  for  the  two  cases  were  similar  over 
smaller  stress  levels,  but  the  GDIS  case  continued  beyond  the  point 
where  the  FSIS  distribution  went  to  zero.  Thus,  the  former  provided 
the  power  law  greater  emphasis  and  damage  rates  were  found,  both 
empirically  and  analytically,  to  be  larger.  Finally,  the  theoreti¬ 
cal  damage  rate  for  the  FSIS  condition  when  Qv  =  380  is  less  than 
unity  while  the  experimental  value  is  greater  than  one  (as  it 
should  be) .  This  result  emphasizes  that  analytical  developments 
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are  valid  only  for  low  Qv  products.  However,  if  information  about 
system  reliability  under  loading  for  high  Qv  products  (or  Qv/co  >  l) 
is  desired,  then  as  suggested  earlier  one  can  safely  use  the  gaussian 
distribution  as  a  guideline. 

In  conclusion,  it  has  been  shown,  on  the  basis  of  a  Palmgren- 
Miner  damage  accumulation  law,  that  a  linear  harmonic  oscillator  is 
much  more  prone  to  fatigue  wearout  for  the  same  response  rms  level 
under  impulsive  noise  excitation  than  for  gaussian  forcing.  This 
bears  out  earlier  damage  estimates,  obtained  with  a  dead-zone  device 
impulsive  noise  generator.  One  safely  may  generalize  that  linear 
system  reliability  decreases  for  impulsive  noise  loading  from  that 
expected  under  gaussian  loading  at  the  same  rms  level;  thus,  blind 
usage  of  the  Rayleigh  peak  distribution  will  generally  lead  to 
dangerously  low  estimates  of  system  damage. 
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Fig.  13  Schematic  Diagram  of  Impulsive  Noise  Generator 


Appendix  I .  The  Impulsive  Noise  Generator 

In  order  to  simulate  impulsive  noise  for  laboratory  study  of 
system  response  it  was  necessary  to  design  a  generator  whose  essen¬ 
tial  characteristics  approximate  those  of  the  theoretical  Poisson 
impulsive  noise  model.  Thus,  some  association  between  the  analyti¬ 
cal  quantities  v  (impulse  frequency)  and  p(a)  (impulse  strength 
distribution)  with  physical  parameters  of  the  generator  must  be 
established . 

Since  it  is  convenient  to  design  switching  circuitry  which 
produces  pulses  whose  widths  are  fixed,  one  may  associate  with  the 
area  of  the  (sufficiently  narrow)  pulse  the  impulse  strength  of 
the  analytical  impulse.  In  this  light,  the  heights  of  the  gener¬ 
ated  pulses  are  solely  indicative  of  the  pulse  strengths.  Thus, 
if  these  pulse  heights  are  modulated  by  a  signal  whose  first-order 
probability  density  is  p(a),  then  the  areas (and  hence  the  impulse 
strengths)  of  the  generated  pulses  are  also  distributed  according 
to  p  (a)  . 


There  still  remains  the  problem  of  triggering  the  pulses  so 
their  occurrences  (in  time)  are  mutually  independent  and  their 
number  Poisson-distributed.  Since  level  crossings  of  a  w^e-band 
gaussian  process  are  believed  to  meet  these  requirements,  °  detec¬ 
tion  by  a  Schmitt  trigger  should  accomplish  this  goal.  If,  in 
turn,  this  detection  circuit  drives  the  aforementioned  switching 
circuitry,  the  generator  is  complete. 

Figure  13  provides  the  detail  for  such  a  generator.  The 
capacitor  C  in  the  monostable  multivibrator  controls  the  pulse- 
duration.  Duration  times  as  short  as  7.3  fj. sec  (a  lower  limit 
determined  by  circuit  rise  time  -  lftsec)  can  be  achieved  and,  for 
all  practical  purposes,  there  is  no  upper  bound  on  the  durations. 
The  impulse  rates  depend  upon  the  frequency  of  occurrence  of  the 
level-crossings  of  the  gaussian  (or  other)  reference  noise  source. 
For  a  pulse-duration  of  500^t.sec  (a  value  used  for  experimental 
studies)  rates  of  from  1  per  second  to  1,000  per  second  were  con¬ 
veniently  obtained. 

Several  checks  were  made  to  determine  how  closely  the  Poisson 
model  is  approximated.  Two  of  them  depend  on  the  ability  of  Camp¬ 
bell’s  Theorem  to  predict  root-mean-square  of  the  lDOF2system 
response.  Equation  (24)  shows  the  response  variance  a  =  varies 
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2  •  •  14 

as  v  and  <a  ).  These  dependences  have  been  investigated  and 

verified  experimentally.  Moreover,  the  time  intervals  between 

level  crossings  of  the  gaussian  source  were  recorded  (up  to  700 

intervals)  and  the  mean  and  variance  computed.  Since  for  an 

exponentially— distributed  time  variable  (e.g.,  time  duration  between 

adjacent  Poisson  process  pulses)  the  standard  deviation  equals  the 

mean,  and  here  also  1/v ,  these  experimentally  measured  quantities 

were  compared  and  found  nearly  identical. 

The  probability  density  of  the  modulator  signal  is  chosen  to 
provide  the  desired  p(a)  in  any  given  situation.  Thus,  for  the 
FSIS  condition  a  zero-mean  square  wave  was  used  while  simulation 
of  the  GDIS  case  required  a  gaussian  input. 
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Appendix  II.  The  Rise  and  Positive  Maxima  Distribution  Analyzer  (RPMDA) 

The  RPMDA  provides  instrumentation  for  the  empirical  measure¬ 
ment  of  rise  (or  fall)  distributions  and  positive  maxima  (or 
negative  minima)  exceedance  probabilities.  In  order  to  obtain  data 
efficiently  both  operations  may  be  performed  simultaneously.  More¬ 
over,  since  both  operations  require  knowledge  of  the  occurrence  of 
a  maximum  (or  minimum) ,  the  detection  equipment  necessary  for  this 
phase  can  be  shared  by  both  devices,  thus  minimizing  the  overall 
complexity  of  the  two  systems. 

The  operation  of  the  RPMDA,  Fig.  14,  can  be  described  by 
following  its  operation  over  one  "cycle",  i.e.,  from  say  one  minimum 
to  the  next.  If  the  input  signal  is  differentiated  and  then  ampli¬ 
fied  by  a  high  gain  amplifier  operating  in  a  saturated  mode,  then 
the  resultant  signal  is  a  constant-magnitude  voltage,  positive  dur¬ 
ing  the  first  "half"  cycle  (the  "half-point"  being  the  time  occur¬ 
rence  of  the  maximum  in  the  cycle)  and  negative  over  the  remaining 
"half".  Thus,  a  square-wave  gate  is  obtained  and  is  shared  by 
both  phases  of  operation  of  the  RPMDA. 

For  rise  measurements  the  above  gate  signal  controls  a  gated- 
integrator  whose  input  is  the  derivative  of  the  process  being 
analyzed.  Hence,  the  output  of  this  integrator  is  the  rise  portion 
of  the  original  process  (within  a  scale  factor)  with  the  minimum 
clamped  to  zero.  An  electronic  counter  placed  at  some  positive 
level  provides  a  detector  for  the  exceedance  of  that  level  and  thus 
the  exceedance  of  the  same  positive  rise  level.  By  changing  the 
detector  level  -  or,  as  is  actually  done,  rescaling  the  input  signal 
to  the  gated  integrator  -  the  number  of  positive  rises  which  exceed 
selected  levels  can  be  counted,  and  the  operation  is  complete. 

Measurement  of  the  maxima  distribution  also  utilizes  the  gate 
signal  discussed  earlier,  but  rather  than  operate  a  gated-integrator 
as  for  the  rise  statistics,  it  is  used  in  a  parallel  fashion  to 
drive  a  gated-amplifier .  Thus,  the  resultant  signal  is  proportional 
to  the  original  process  (note,  here  the  minimum  is  not  clamped  to 
zero,  but  is  at  its  true  level  of  occurrence)  during  the  first  "half" 
of  the  cycle  and  zero  otherwise .  .\nother  electronic  counter  at  a 
positive  level  provides  the  detection  scheme  as  before,  and  as  the 
scaling  of  the  input  to  the  gated-amplifier  is  changed,  the  number 
of  positive  maxima  occurring  above  that  level  as  a  function  of  that 
level  is  obtained,  i.e.,  the  exceedance  probability  of  positive 
maxima. 
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Fig.  14  Schematic  Diagram  of  the  Rise  and  Positive  Maxima  Distribution  Analyzer 


Figure  14  shows  a  schematic  diagram  utilizing  four  operational 
amplifiers  (two  Type-0  Tektronix  plug-in  units) .  Since  wide-band 
differentiation  is  difficult,  if  not  impossible,  the  device  is 
designed  to  function  on  band-limited  processes  with  the  ratio  of 
lower  to  upper  cutoff  frequency,  equal  to  1/20.  This  restric¬ 

tion  is  coupled  via  equipment  limitations  to  the  dynamic  range  of 
measurements  allowed;  for  the  equipment  used,  rises  and  maxima 
between  0.2  and  10  times  rms  can  be  detected.  The  allowed  frequency 
band  can  be  changed  simply  by  altering  the  element  C  ,  C  ,  R  , 
and  C  .  D  N  E 

A  detailed  discussion  of  modifications  imposed  on  the  afore¬ 
mentioned  concepts  by  design  considerations,  and  also  the  calibra¬ 
tion  procedures  and  parametric  values  for  various  bandwidth  of  the 
device,  is  available  in  an  internal  report  by  M.  Joselevich  and 
G .  L .  Hedin . 
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